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ABSTRACT 

Global VLBI (EVN+VLBA) polarization observations at 5 and 8.4 GHz of ten high 
redshift {z > 3) quasars are presented. The core and jet brightness temperatures 
are found through modelling the self-calibrated Mu-data with Gaussian components, 
■ which provide reliable estimates of the flux density and size of individual components. 

The observed high core brightness temperatures (median Tb,corc = 4 x 10^^ K) are 
consistent with Doppler boosted emission from a relativistic jet orientated close to 
the line-of-sight. This can also explain the dramatic jet bends observed for some of 
our sources since small intrinsic bends can be significantly amplified due to projection 
effects in a highly beamed relativistic jet. We also model-fit the polarized emission 
and, by taking the minimum angle separation between the model-fitted polarization 
angles at 5 and 8.4 GHz, we calculate the minimum inferred Faraday rotation mea- 
sure (RMniin) for each component. We also calculate the minimum intrinsic RM in 
the rest frame of the AGN (RMj^?^ = RMi„in(l + zf), first subtracting the integrated 
(presumed foreground) RM in those cases where we felt we could do this reliably. The 
resulting mean core |RM^;^| is 5580 rad m~^, with a standard deviation of 3390 rad 
m~^, for four high-z quasars for which we believe we could reliably remove the fore- 
ground RM. We find relatively steep core and jet spectral index values, with a median 
core spectral index of —0.3 and a median jet spectral index of —1.0. Comparing our 
results with RM observations of more nearby Active Galactic Nuclei at similar emit- 
ted frequencies does not provide any significant evidence for dependence of the quasar 
nuclear environment with redshift. However, more accurate RM and spectral informa- 
tion for a larger sample of sources would be required before making any conclusive 
statements about the environment of quasar jets in the early universe. 

Key words: radio continuum: galaxies - galaxies: jets - galaxies: magnetic fields 



1 INTRODUCTION 

Quasars at 2: > 3 and with GHz flux densities of the order of 
1 Jy have luminosities of ~ 10^* W Hz~^, making them the 
most luminous, steady radio emitters in the Universe. In this 
paper, we present 5 and 8.4 GHz global VLBI polarization 
observations of ten quasars at 2 > 3. There have been sev- 
eral recent studies of high-redshift je t s on VLBI scales (e.g., 
iGurvits et alll2000l : IfVev et all [200^ : I Veres et al.ll2010l ') but 
our observations are the first to explore this region of "lumi- 
nosity - emitted frequency" parameter space in detail with 
polarization sensitivity. In the standard ACDM cosmology 
(Ho = 70 km s"^ Mpc"\ = 0.7 and VLm = 0.3), the an- 
gular diam eter distance reaches a maximum at z ~ 1.6 (e.g., 
lHogel[l999h . meaning that more distant objects actually be- 



gin to increase in apparent angular size with increasing 2. 
As a consequence, the linear scale for observations of ob- 
jects at 2 ~ 3 is similar to objects at 2 ~ 0.7. Therefore, our 
measurements allow us to study pc-scale structures at emit- 
ted frequencies of ~20-45 GHz in comparison with prop- 
erties of low-redshift core-dom inated Active Galactic Nucle i 
(AGN) at 22 and 43 GHz (e.g.. lO'SulIivan fc Gabu^l2009l ) 
at matching length scale and emitted frequency. 

An important question we try to address with these 
observations is whether or not quasar jets and their sur- 
rounding environments evolve with cosmic time. Any sys- 
tematic diiferences in the polarization properties or the fre- 
quency dependence of the polarization between high and low 
redshift sources would suggest that conditions in the cen- 
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tral engines, jets, and/or surrounding media of quasars have 
evolved with redshift. Previous observations of h igh redshift 
quasars fe.g. jFrev et al.l[l997l : IParagi et al.lll999f ) often show 
that these quasars can be strongly dominated by compact 
cores and that their extremely high apparent luminosities 
are likely due to large Doppler boosting, implying that their 
jets move at relativistic velocities and are pointed close to 
our line of sight. With this geometry, it is difficult to distin- 
guish emerging jet components from the bright core. How- 
ever, high resolution polarization observations have proven 
very effective in identifying barely resolved but highly po- 
larized jet components, whose intensity is swamped by the 
intensity of the nearby core, but whose polarized flux is 
comparable to or even greater than that of the core (e.g., 
iGabuzda 19991 . Thus, polarization sensitivity can be crucial 
to understanding the nature of compact jets in the highest 
redshift quasars. 

There is also some evidence that the VLBI jets of high 
redshift quasars are commonly bent or distorted; one re- 
mark able case revealed b y VSOP observations is of 1351- 
018 l|Gurvits et al.l I2OO0I I where the jet appears to bend 
through almost 180°. Polarization observations can help 
elucidate the physical origin of these observed changes in 
the jet direction in the plane of the sky since the direc- 
tion of the jet polarization can refl ect the direction of the 
underlying flow. In some q uasars jCawthorne et al.l 1 19931 : 
lO'Sullivan fc Gabuzda|[2009l ) , the observed polarization vec- 
tors follow curves in the jet, indicating that the curves are 
actual physical bends, while in others, remarkable unifor- 
mity of polarization position angle has been observed over 
curved sections of jet, suggesting that the apparent curves 
correspond to bright regions in a much broader underly- 
ing flow. One possible origin for bends could be collisions 
of the jet with denser areas in the surrounding medium as 
suggested, for ex ample, in the case o f the jet in the radio 
galaxy 4C 41. 17 jGurvits et al.lll997] 'l. The polarization in- 
formation allows us to test for evidence of such interaction, 
where the Faraday rotation measure (RM) and degree of 
polarization should increase substantially in those regions. 

In Section 2, we describe our observations and the data 
reduction process. Descriptions of the results for each source 
are given in Section 3, while their implications for the jet 
structure and environment are discussed in Section 4. Our 
conclusions are presented in Section 5. We use the S oc v" 
definition for the spectral index. 



2 OBSERVATIONS AND DATA REDUCTION 

Ten high-redshift {z > 3) AGN jets (fisted in Table[l]), which 
were all previously successfully imaged using global VLBI 
baselines, were targeted for global VLBI polarization obser- 
vations at 4.99 and 8.415 GHz with two IPs of 8 MHz band- 
width at each frequency. These sources do not comprise a 
complete sample in any sense, and were chose n base d on the 
previous VLBI observations of iGurvits et al.l (|200Cll ). which 
indicated the presence of components with sufficiently high 
total intensities to suggest that it was feasible to detect their 
polarization. 

The sources were observed for 24 hours on 2001 June 5 
with all 10 antennas of the NRAO Very Long Baseline Array 
(VLBA) plus six antennas of the European VLBI Network 



Table 2. List of telescopes included in the global VLBI experi- 
ment. 



Telescope Location Code Diameter Comment 

[m] 



EVN 



Effelsberg, Germany 


EB 


100 




Medicina, Italy 


MC 


32 




Noto, Italy 


NT 


32 




Torun, Poland 


TR 


32 


Failed 


Jodrell Bank, UK 


JB 


25 




Westerbork, Netherlands 


WB 


93* 




VLBA, USA 


Saint Croix, VI 


SC 


25 




Hancock, NH 


HN 


25 




North Liberty, lA 


NL 


25 




Fort Davis, TX 


FD 


25 




Los Alamos, NM 


LA 


25 


Failed 


Pietown, NM 


PT 


25 




Kitt Peak, AZ 


KP 


25 




Owens Valley, CA 


OV 


25 




Brewster, WA 


BR 


25 




Mauna Kea, HI 


MK 


25 




* Equivalent diameter for 


the phased 


array, y/lA 


X 25 m. 



(EVN) capable of fast frequency switching between 5 and 
8.4 GHz (see Table [2] for a list of the telescopes used). The 
target sources were observed at alternating several-minute 
scans at the two frequencies. This led to practically simulta- 
neous observations at the two frequencies with a resolution 
of ~ 1.1 mas at 5 GHz and ~ 0.7 mas at 8.4 GHz (for the 
maximum baseline of 11,200 km from MK to NT). 

The data were recorded at each telescope with an aggre- 
gate bit rate of 128 Mbits/s, recorded in 8 base-band chan- 
nels at 16 Msamples/sec with 2-bit sampling, and correlated 
at the Joint Institute for VLBI in Europe, Dwingeloo, The 
Netherlands. Los Alamos (LA) was unable to take data due 
to communication problems and Toruh (TR) also failed to 
take data, so both were excluded from the processing. 

The 5 and 8.4 GHz VLBI data were calibrated inde- 
pendently using standard techniques in the NRAO AIPS 
package. In both cases, the reference antenna for the VLBI 
observations was the Kitt Peak telescope. The instrumen- 
tal polarizations ('D-terms') at each frequency were deter- 
mined using observations of 3C 84, using the task LPCAL 
and assuming 3C 84 to be unpolarized. The polarization an- 
gle (2x = arctan^) calibration was achieved by comparing 
the total VLBI-scale polarizations observed for the compact 
sources 1823+568 (at 5 GHz) and 2048-f312 (at 8.4 GHz) 
with their polarizations measured simultaneously with the 
NRAO Very Large Array (VLA) at both 5 and 8.4 GHz, 
and rotating the VLBI polarization angles to agree with the 
VLA values. We obtained 2 hours of VLA data, overlapping 
with the VLBI data on June 6. The VLA angles were cali- 
brated using the known polarization angle of 3C 286 at both 
5 and 8.4 GHz (see AIPS cookbool^, chapter 4). At 5 GHz, 



http: / /www. aips.nrao.edu/cook. html 
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Table 1. Basic information for the high-rcdshift quasars and VLBI data presented in this paper. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Name 


z 


Linear Scale 


Frcq. 


Beam 


BPA 


Aotal 


^poak 


<^/,rms 








[pc mas ^] 


[GHz] 


[mas X mas] 


[dog] 


[mJy] 


[mjy/bm] 


[mJy/bm] 


[10^* W Hz"!] 


0014+813 


3.38 


7.55 


4.99 


1.10 X 0.95 


-85.9 


861.4 


583.2 


0.3 


2.10 








8.41 


0.95 X 0.68 


-85.1 


645.1 


445.6 


0.9 


1.57 


0636+680 


3.17 


7.71 


4.99 


1.25 X 0.91 


-1.3 


348.8 


270.0 


0.3 


0.76 








8.41 


0.82 X 0.62 


5.3 


281.4 


162.7 


0.5 


0.62 


0642+449 


3.41 


7.53 


4.99 


1.70 X 0.85 


-15.6 


2248.2 


1974.9 


0.7 


5.55 








8.41 


1.03 X 0.70 


14.4 


2743.9 


2460.4 


3.9 


6.77 


1351-018 


3.71 


7.30 


4.99 


3.51 X 1.10 


-7.0 


788.0 


662.2 


0.4 


2.23 








8.41 


2.07 X 0.86 


-2.8 


594.0 


534.2 


1.4 


1.68 


1402+044 


3.21 


7.68 


4.99 


3.52 X 0.95 


-7.7 


831.7 


617.4 


0.5 


1.86 








8.41 


2.03 X 0.75 


-6.8 


719.7 


488.2 


0.3 


1.61 


1508+572 


4.30 


6.88 


4.99 


1.23 X 0.97 


15.2 


291.1 


248.7 


0.2 


1.04 








8.41 


0.80 X 0.76 


8.9 


198.3 


170.3 


0.5 


0.71 


1557+032 


3.90 


7.16 


4.99 


3.81 X 0.97 


-8.2 


287.0 


249.6 


0.3 


0.88 








8.41 


2.06 X 0.80 


-8.8 


255.8 


239.1 


1.0 


0.78 


1614+051 


3.21 


7.68 


4.99 


3.74 X 0.98 


-8.1 


792.9 


588.9 


0.5 


i.rr 








8.41 


1.99 X 0.80 


-9.4 


465.0 


333.6 


1.0 


1.04 


2048+312 


3.18 


7.70 


4.99 


2.17 X 1.12 


-21.2 


547.0 


192.7 


1.5 


1.20 








8.41 


2.06 X 0.71 


-16.2 


539.1 


253.9 


0.6 


1.19 


2215+020 


3.55 


7.42 


4.99 


3.45 X 1.03 


-7.6 


271.0 


190.6 


0.2 


0.71 








8.41 


2.06 X 0.71 


-6.8 


174.3 


122.6 


0.3 


0.46 



Column designation: 1 - source name (lAU B1950.0); 2 - rcdshift; 3 - projected distance, in parsecs, corresponding to 1 mas on the sky; 4 - observing 
frequency in GHz; 5 - size of convolving beam in mas X mas at the corresponding frequency for each source; 6 - beam position angle (BPA), in degrees; 
7 - total flux density, in mJy, from CLEAN map; 8 - peak brightness, in mJy beam" ; 9 - brightness rms noise, in mJy beam" , from off-source region; 
10 - monochromatic total luminosity at the corresponding frequency in units of 10^^ W Hz""*". 



we found A^sghz = Xvla — xvlbi = 67.5° and at 8.4 GHz, 
Ax8.4GHz = Xvla — Xvlbi = 83.5°. 



ponent size reflects the t rue size o f that particular jet emis- 
sion region or not. From iLobanovl (|2005D . we have 



3 RESULTS 

Table [1] lists the 10 high-z quasars, their redshifts, the 
FWHM beam sizes at 5 and 8.4 GHz, the total CLEAN 
flux and the peak flux densities of the 5 and 8.4 GHz maps, 
the noise levels in the maps, and the 5 and 8.4 GHz lumi- 
nosities, calculated assuming a spectral index of zero and 
isotropic radiation. Since the emission from these sources is 
likely relativistically boosted, the true l uminosity is proba - 
bly smaller by a factor of 10-100 (e.g., ICohen et aLl^2007^ . 
Figure [T] gives an example of the sparse but relatively uni- 
form uv coverage obtained for this experiment. The 5-GHz 
uv coverage for 2048-1-312 is shown; the 8.4 GHz is essen- 
tially identical but scaled accordingly. 

The DIFMAP software package (|Shepherdl 119971 ) was 
used to fit circular and/or elliptical Gaussian components 
to model the self-calibrated source structure. The bright- 
ness temperatures in the source frame were calculated for 
each component from the integrated flux and angular size 
according to 



Tt = 1.22 X 10 



12 >S'total(l + z) 



dimLxdm.ini' 



(1) 



where the total fiux density is measured in Jy, the FWHM 
size is measured in mas and the observing frequency is mea- 
sured in GHz. The limiting angular size for a Gaussian com- 
ponent (diim) was also calculated to check whether the com- 



, /ln2, , , f S/N 

"lim — \l — ^OmaxOmin m I 



(2) 



where 6max and femin are the major and minor axes of the 
FWHM beam and S/N is the signal-to- noise ratio of a par- 
ticular component. In other words, component sizes smaller 
than this value yielded by the model-fitting are not reliable. 
Table[3]shows the results of the total intensity model-fltting 
for each source. For clarity of comparison, the component 
identified as the core has been defined to be at the origin, 
and the positions of the jet components determined rela- 
tive to this position. The errors in the model-fitted posi- 
tions were estimated as Ar = Urmsd-^l + S'poak/o-rms/2S'poak 
and A9 — arctan(Ar/r), where arms is the residual noise of 
the map after the subtraction of the model, d is the model- 
fitted componen t size and Speak is the peak flux density (e.g., 
[Pomalont 1999i: iLee et~al . 2008; Kudryavtscva ct al. 20iq). 
These are formal errors, and may yield position errors that 
are appreciably underestimated in some cases; when this oc- 
curs we use a la error estimate of ±0.1 beamwidths in the 
structural position angle of the component in question. 

Polarization model flts, listed in Table 4, were found 
using the Brandeis VISFIT package (|Roberts et al.lll987l : 
iGabuzda et al.l 19891) adapted to run in a linux environment 
bv iBezrukovs fc Gabuzdal (|2006l '). The positions in Table 4 
have been shifted in accordance with our cross-identification 
of the corresponding intensity cores, when we consider this 
cross-identification to be reliable. The errors quoted for the 
model-fitted components are formal la errors, correspond- 
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ing to an increase in the best fitting by unity; again, we 
have adopted a la error estimate of ±0.1 beamwidths in the 
structural position angle of the component in question when 
the position errors are clearly underestimated. 

Before constructing spectral-index maps, images with 
matched resolution must be constructed at the two frequen- 
cies. Since the intrinsic resolutions of the 8.4 and 5-GHz 
images were not very different (less than a factor of two), 
we achieved this by making a version of the final 8.4-GHz 
image with the same cell size, image size and beam as the 
5-GHz image. The two images must also be ahgned based on 
optically thin regions of the structure: the mapping proce- 
dure aligns the partially optically thick cores with the map 
origin, whereas we expect shifts between the physical posi- 
tions of the cores at the two frequencies. When possible, we 
aligned the two images by comparing the positions of opti- 
cally thin jet components at the two frequencies derived from 
model-fits to the VLBI data, or using the cross-correlation 
technique of Croke & Gabuzda (2008). After this alignment, 
we constructed spectral-index maps in AIPS using the task 
COMB. 

In the absence of Faraday rotation, we expect the polar- 
ization angles for corresponding regions at the two frequen- 
cies to coincide; in the presence of Faraday rotation of the 
electric vector position angle (EVPA), which occurs when 
the polarized radiation passes through regions of magnetized 
plasma, the observed polarization angles will be subject to a 
rotation by RMA^, where RM is the Faraday rotation mea- 
sure and A is the observing wavelength. We are not able to 
unambiguously identify the action of Faraday rotation based 
on observations at only two frequencies; however, Faraday 
rotation provides a simple explanation for differences in the 
polarization angles observed at the two frequencies. We ac- 
cordingly calculated tentative Faraday rotation measures, 

RMmi„ = (X5 GHz - X8.4 GHz)/(Ai GHz " ^8.4 GHz), based 

on the results of the polarization model fitting, taking the 
minimum difference between the component EVPAs. We 
also constructed tentative RM maps in AIPS using the task 
COMB. In both cases, if Faraday rotation is operating, the 
derived RM values essentially represent a lower limit for the 
true Faraday rotation (assuming an absence of nvr rotations 
in the observed angles). 

With our two frequencies, the nvr ambiguity in the RM 
corresponds to n times 1350 rad m~^. Table 4 lists the values 
of RMmin and the corresponding intrinsic polarization angle 
Xo, obtained by extrapolating to zero wavelength. However, 
we emphasize that polarization measurements at three or 
more wavelengths are required to come to any firm conclu- 
sions about the correct RM and xo values for these sources. 

The main source of uncertainty for the RMs (apart from 
possible nvr ambiguities) comes from the EVPA calibration, 
which we estimate is accurate to within ±3°; this corre- 
sponds to an RM error of ±23 rad vcC^ between 5 and 
8.4 GHz. We have also attempted to obtain the correct 
sign and magnitude of the RM in the immediate vicinity 
of the AGNs by subtracting the integrated RMs derived 
from lower-frequency VLA measurements centered on our 
sources, obtained from the literature. Note that these rep- 
resent the integrated RMs directly along the line of sight 
toward our sources, rather than along a nearby sight-line; 
the typical uncertainties in such measurements when based 
on observations at several frequencies near 1-2 GHz are typ- 



Table 5. Total Intensity and Polarization image cutoff values. 



5 GHz 8.4 GHz 



Source name 


/ cutoff 


p cutoff 


I cutoff 


p cutoff 




[mjy/bm] 


[mJy/bm] 


[mJy/bm] 


[mJy/bm] 


0014+813 


1.2 


0.7 


2.8 


1.3 


0636+680 


0.9 


0.7 


1.5 


0.2 


0642+449 


3.6 


2.4 


5.0 


11.0 


1351-018 


1.3 


1.1 


4.0 


1.3 


1402+044 


2.0 


1.4 


1.1 


1.7 


1508+572 


1.1 


0.7 


1.5 


1.8 


1557+032 


1.1 


0.9 


1.0 


1.2 


1614+051 


1.7 


1.1 


4.2 


1.5 


2048+312 


6.0 


1.9 


1.9 


1.5 


2215+020 


0.7 


0.7 


1.6 


0.9 



ically no larger than about ±5 rad/m^. Because of the lower 
resolution of these measurements and the greater promi- 
nence of the jets at lower frequencies, the polarization de- 
tected in such observations usually originates fairly far from 
the VLBI core, where we expect the overall RM local to the 
source to be negligible. Therefore, the lower-frequency inte- 
grated RMs usually correspond to the foreground (Galactic) 
contribution to the overall RM detected in our VLBI data, 
and subtracting off this value should typically help isolate 
the RM occurring in the immediate vicinity of the AGN. 
Due to the extremely high redshifts of these sources it can 
be very important to remove the foreground RM before esti- 
mating the intrinsic minimum RM in the source rest frame, 
RMiSf^ = (1± 2)2 (RMmin -RMcai). Table 4 also lists RMcai 
and RMJ^i^ for sources which have measured integrated RMs 
in the literature. We take the polarization angles measured 
at the two frequencies to essentially be equal (negligible RM) 
within the errors if they agree to within 6°, and do not at- 
tempt to correct such angles for Galactic Faraday rotation, 
since there can be some uncertainty about interpreting in- 
tegrated RMs as purely Galactic foreground RMs, and the 
subtraction of the integrated values from such small nomi- 
nal VLBI RM values could lead to erroneous results in some 
cases. 

For each source, we present the total-intensity distribu- 
tions at 5 and 8.4 GHz overlaid with EVPA sticks in Fig. 2 
and also plot the positions of the model-fitted total intensity 
and polarization components. The derived spectral-index 
and RMmin distributions generally do not show unexpected 
features, and are not presented here; these can be obtained 
by contacting S.P. O' Sullivan directly. 

3.1 0014+813 

This source was discovered in the radio and cla s sified as a flat 
spectr um radio quasar bv iKuehr et al.1 (|l98ll ). iKuehr et al.l 
(|l983t ) obtained a redshift of 3.366, and found the source 
to be exceptionaUy brig ht in the optical but unpolarized, 
while iKaspi et al.1 (|2007l ) found very little optical variability 
for th e source. From VL BA observations at 8 GHz over 5 
years, IPiner et al.1 |20o3) did not detect any outward mo- 
tion. The source was observed by the Swift satellite in 2007 
January from optical to hard X-rays. Thro ugh construction 
of the spectral energy distribution (SED), IChisellini et al.l 
(2009) show that it may harbour one of the largest black hole 
masses ever inferred. By associating the strong optical-UV 
emission with a thermal origin from a standard optically- 



© 2010 RAS, MNRAS 000, [T]-?? 



VLBI Polarization of High-z Jets 5 



Table 3. Total intensity model fit parameters of all sources. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Name 


Frcq. 


Comp. 


r 


9 




dmax 


d. 






Tb 




[GHz] 




[mas] 


[dog] 


[m.Jy] 


[mas] 


[mas] 


[mas] 


[XIO"' K] 




4.99 










0.53 





.32 





,03 


78.2 + 6.4 






B 


0.66 ± 0.10 


— 175.7 ± 1.2 


196 ± 11 


0.53 





.53 





.05 


15.2 + 2.3 






C 


4.93 ± 0.10 


-162.8 ± 1.1 


13 ± 3 


0.92 





.92 





,21 


0.32 + 0.21 






D 


6.41 ± 0.21 


-171.0 ± 1.9 


16 ± 4 


1.75 


1 


,75 


0, 


,19 


0.11 + 0.08 






E 


9.51 ± 0.16 


-170.9 ± 1.0 


17 ± 4 


1.49 


1 


,49 


0, 


,18 


0.16 + 0.11 




8.41 


A 






547 ± 36 


0.46 





13 





04 


65.9 + 11.0 






B 


0.74 ± 0.07 


-178.3 ± 1.7 


85 ± 14 


0.40 


0, 


,40 


0, 


,09 


3.9 + 1.7 


0636+680 


4.99 


A 






344 ± 15 


0.69 





35 





03 


28.9 + 3.1 






B 


2.22 ± 0.54 


140.3 ± 9.0 


6 ± 3 


2.76 


2 


.76 





.23 


0.02 + 0.02 




8.41 








ZOO ± 21 


0. 78 


0, 


.38 


0, 


.04 


7.4+1.4 


0642+449 


4.99 


A 






2209 ± 60 


0.41 


0, 


,28 


0, 


,02 


417 + 28 






B 


3.22 ± 0.45 


92.5 ± 7.5 


45 ± 14 


2.86 


2 


.86 


0, 


.15 


0.12 + 0.11 




8.41 


A 


- 


- 


2313 ± 134 


0.16 





.16 





.03 


731 ± 103 






Al 


0.20 ± 0.09 


158.1 ± 0.8 


458 ± 60 


0.09 


0, 


.09 


0, 


.08 


439 ± 138 


1351-018 


4.99 


A 


- 


- 


681 ± 23 


0.77 





.23 





.05 


89.5 + 7.5 






B 


1.01 ± 0.16 


132.7 ± 1.0 


79 ± 8 


0.48 





,48 





,14 


7.9 + 2.0 






C 


1.64 ± 0.12 


56.6 ± 0.5 


14 ± 3 


0.17 


0, 


,17 


0, 


,32 


- 






D 


11.85 ± 0.59 


-12.9 ± 2.8 


17 ± 5 


4.02 


4, 


,02 


0, 


,30 


0.02 + 0.02 




8.41 


A 






578 ± 37 


0.67 


0, 


,06 


0, 


,06 


125.9 + 20.1 






B 


0.83 ± 0.12 


138.5 ± 9.7 


13 ± 6 


0.76 


0, 


,76 


0, 


,40 


0.18 + 0.19 






C 


5.96 ± 0.63 


-2.1 ± 5.8 


3 ± 3 


2.43 


2 


,43 


1 


,06 


0.003 + 0.008 


1402+044 


4.99 


A 


- 


- 


330 ± 17 


0.23 


0, 


,23 


0, 


,05 


130.0 + 16.7 






B 


0.76 ± 0.25 


-24.5 ± 1.3 


374 ± 19 


0.46 


0, 


,46 


0, 


,04 


36.9 + 4.6 






C 


3.33 ± 0.15 


-44.4 ± 4.5 


19 ± 5 


1.93 


1 


,93 


0, 


,23 


0.11 + 0.08 






D 


9.07 ± 0.14 


-46.6 ± 0.9 


64 ± 9 


2.13 


2 


,13 





,11 


0.29 + 0.12 








12. D4 ± 0. 00 


— II./ ± 2..3 


44± 10 




4 


,72 





,18 


0.04 It O.Oo 




8.41 


A 






300 ± 19 


0.14 





,14 





,07 


111.1 + 17.5 








O.oo ± 0. ID 


— 2d. ^ It i .U 


03D It 21 


0.43 





,43 





06 








C 


3.23 ± 0.29 


— 43.3 ± 5.5 


13 ± 5 


1.62 


1 


,62 


0, 


,27 


0.03 ± 0.04 






D 


9.22 ± 0.20 


-46.2 ± 1.2 


55 ± 11 


2.07 


2 


,07 


0, 


,15 


0.09 ± 0.05 






E 


12.83 ± 0.32 


-77.8 ± 1.5 


20 ± 7 


2.09 


2 


,09 


0, 


18 


0.03 + 0.03 


1508+572 


4.99 


A 


_ 


_ 


276 ± 11 


0.45 


0, 


,28 


0, 


,03 


56.1 + 5.3 






B 


2.00 ± 0.12 


-178.4 ± 3.3 


14 ± 3 


1.37 


1 


,37 


0, 


13 


0.20 + 0.11 




8.41 


A 






196 ± 14 


0.42 


0, 


14 


0, 


,03 


30.0 + 5.2 






B 


1.94 ± 0.14 


-172.1 ± 3.7 


6 ± 3 


0.70 


0, 


,70 


0, 


,20 


0.11 + 0.13 


1557+032 


4.99 


A 






262 ± 12 


0.70 





,26 





,06 


34.0 + 3.9 






B 


2.64 ± 0.17 


139.4 ± 2.4 


14 ± 3 


1.24 


1 


24 





26 


0.22 + 0.13 




8.41 


A 






252 ± 15 


0.49 


0, 


12 


0, 


,05 


35.0 + 5.1 






B 


2.16 ± 0.37 


148.5 ± 9.2 


15 + 5 


2.46 


2 


,46 





,18 


0.02 + 0.02 


1614+051 


4.99 


A 






566 ± 24 


0.37 


0, 


,37 


0, 


,05 


63.5 + 6.7 






B 


1.26 ± 0.16 


— 154.4 ±0.7 


229 ± 15 


0.54 


0, 


54 


0, 


,08 


12.4 + 2.1 




8.41 


A 


- 


- 


347 ± 28 


0.24 


0, 


,24 


0, 


,07 


32.4 + 6.5 






B 


1.26 ± 0.12 


— 154.1 ± 0.7 


117 ± 17 


0.37 


0, 


,37 


0, 


12 


4.6 +1.6 


2048+312 


4.99 


A 


- 


- 


393 ± 38 


1.78 


1 


,24 


0, 


,10 


3.7 + 1.0 






B 


1.55 ± 0.11 


65.7 ± 2.2 


126 ± 21 


1.43 


1 


,43 


0, 


,17 


1.3 + 0.6 






C 


8.86 ± 0.36 


56.8 ± 2.0 


20 ± 10 


1.68 


1 


,68 


0, 


44 


0.14 + 0.19 




8.41 


A 






411 ± 23 


0.80 


0, 


,44 


0, 


,04 


8.4 + 1.2 






B 


1.49 ± 0.07 


80.6 ± 2.6 


129 ± 14 


1.29 


1 


,29 


0, 


,07 


0.56 + 0.18 


2215+020 


4.99 


A 






141 ± 9 


0.26 


0, 


,26 


0, 


,08 


45.9 + 7.3 






B 


0.52 ± 0.11 


108.0 ± 2.2 


73 + 7 


0.39 


0, 


,39 


0, 


12 


10.7 + 2.4 






C 


1.53 ± 0.10 


79.1 ± 5.4 


11+3 


1.26 


1 


,26 


0, 


,30 


0.15 + 0.10 






D 


50.52 ± 1.52 


79.1 ± 1.6 


33 ± 11 


9.17 


9, 


.17 


0, 


.17 


0.009 + 0.009 






E 


57.71 ± 2.48 


74.7 ± 2.3 


19 + 8 


11.26 


11.26 


0, 


23 


0.003 + 0.004 




8.41 


A 






94 + 7 


0.20 


0.20 


0, 


,06 


19.0 + 3.5 






B 


0.50 ± 0.08 


115.7 ± 1.3 


66 + 6 


0.22 


0.22 


0, 


,07 


10.7 + 2.4 






C 


0.95 ± 0.07 


95.5 ± 3.6 


17 + 3 


0.64 


0.64 


0, 


14 


0.32 + 0.2 



Column designation: 1 - source name (lAU B1950.0); 2 - observing frequency, in GHz; 3 - component identification; 4 - distance of component from core, in 
mas; 5 - position angle of component with respect to the core, in degrees; 6 - total flux of model component, in mjy; 7 - FWHM major axis of Gaussian 
component, in mas; 8 - FWHM minor axis of Gaussian component, in mas; 9 - Minimum resolvable size, in mas; 10 - measured brightness temperature in 

units of 10^"k. 
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Table 4. Polarization model fit parameters of all sources. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) (12) 


Name 


Frcq. 


Comp. 


r 


e 


P 


X 


m 




xo 






[GHz] 




[mas] 


[dog] 


[mJy] 


[dcg] 


[%] 


[rad m~'^] 


[dog] 


[rad m~^] [rad m~*^] 


0014+813 


4.99 


A 
B 


4.72 ± 0.10 


-163.5 ± 2.1 


4.5 ± 0.3 
1.3 ± 0.5 


-40.9 
47.7 


1.1 ± 0.3 
20.0 ± 4.5 


384 


-120 ± 13 


9" 7194 ±525 




8.41 


A 






2.3 ± 0.2 


-92.2 


0.7 ± 0.2 








0636+680 


4.99 
8.41 
4.99 


A 






4.1 ± 0.3 


78.4 


1.7 ± 0.5 








0642+449 


A 






21.2 ± 1.8 


15.6 


1.1 ± 0.4 










8.41 


A 






68.2 ± 4.2 


15.3 


2.6 ± 0.7 








1351-018 


4.99 


A 
B 


1.45 ± 0.11 


100.3 ± 2.7 


5.2 ± 0.4 
2.8 ± 0.7 


-76.2 
52.3 


0.9 ± 0.3 
2.7 ± 0.5 


-223 
199 


-30 ± 6 
11 ± 3 


-8° -4770 ± 599 
-8° 4592 ± 646 




8.41 


A 
B 


2.38 ± 0.16 


155.9 ± 2.4 


4.1 ± 0.4 
1.7 ± 0.5 


-46.4 
25.7 


0.5 ± 0.2 
1.3 ± 0.4 








1402+044 


4.99 


A 
B 
C 
D 
B 
F 


0.77 ± 0.11 
1.50 ± 0.12 
8.12 ± 0.15 
9.56 ± 0.14 
13.41 ± 0.10 


— 61.5 ± 1.8 
-60.9 ± 3.3 
-45.4 ± 0.9 
-49.2 ± 1.2 
-81.9 ± 1.7 


3.1 ± 0.3 
11.2 ±1.1 
1.6 ± 0.5 
1.1 ± 0.4 
2.0 ± 0.5 
1.6 ± 0.5 


15.7 
4.1 
-26.2 
83.8 
-1.5 
49.8 


1.8 ± 0.5 

2.0 ±0.7 

2.9 ± 0.8 
17.8 ± 5.4 

9.1 ± 2.1 
21.2 ± 6.2 


-199 
— 57 
141 
-59 

- 


57 ± 12 
16 ± 12 
-55 ± 15 
96 ± 49 

- 


- 
- 




8.41 


Al 
A 
B 

c 

D 


2.77 ± 0.20 

0.54 ± 0.09 
1.48 ± 0.08 
8.26 ± 0.11 


177.5 ± 3.4 

-60.1 ± 2.7 
-72.3 ± 2.4 
-43.0 ± 0.9 


5.8 ± 0.6 

2.3 i 0.5 

7.4 ± 0.7 
2.8 ± 0.6 
2.8 ± 0.6 


-6.2 
42.3 
11.7 
-45.1 
91.7 


3.2 ± 0.7 
1.5 ± 0.3 

1.3 ± 0.3 
1.2 ± 0.3 
22.9 ±9.1 








1508+572 


4.99 


A 
B 


2.36 ± 0.18 


177.3 ± 1.4 


7.5 ± 0.9 
0.6 ± 0.3 


-89.4 
36.1 


2.9 ± 0.8 
16.9 ± 5.7 


- 


- 






8.41 


A 


- 


- 


4.0 ± 1.2 


-93.4 


1.3 ± 0.4 








1557+032 


4.99 


A 
B 


1.77 ± 0.19 


19.5 ± 3.6 


3.4 ± 0.5 
0.7 ± 0.3 


54.2 
85.2 


1.7 ± 0.5 

1.8 ± 0.5 


—49 
68 


64 ± 40 
71 ± 32 


3*^ — 1249 ± 713 
3° 1561 ± 642 




8.41 


A 
B 


1.25 ± 0.09 


59.1 ± 1.7 


3.5 ± 0.4 
2.1 ± 0.3 


60.8 
76.2 


1.2 ± 0.4 
6.0 ± 0.7 








1614+051 


4.99 


C 


1.51 ± 0.10 


-95.7 ± 1.0 


3.8 ± 0.5 


-38.5 


1.8 ± 0.5 


522 


-107 ± 10 


8^^ 9110 ±488 




8.41 


A 
B 
C 


0.82 ± 0.08 
1.55 ± 0.09 


-119.1 ± 2.9 
-121.9 ± 3.2 


2.9 ± 0.6 

6.5 ± 0.9 

2.6 ± 0.5 


62.4 
78.7 
-69.8 


3.4 ± 0.7 
2.1 ± 0.4 
2.0 ± 0.3 








2048+312 


4.99 


A 
B 


2.31 ± 0.11 


88.1 ± 1.6 


3.8 ± 0.6 
1.1 ± 0.5 


-39.5 
-61.2 


3.3 ± 0.2 
4.7 ± 0.4 










8.41 


A 
B 


2.04 ± 0.09 


86.6 ± 1.9 


5.5 ± 0.7 
1.4 ± 0.5 


-35.0 
-93.6 


2.2 ± 0.3 
6.2 ± 0.6 








2215+020 


4.99 


A 






0.6 ± 0.2 


57.0 


0.9 ± 0.4 










8.41 


A 






1.5 ± 0.4 


53.8 


1.3 ± 0.5 









Column designation: 1 - source name {lAU B1950.0); 2 - observing frequency, in GHz; 3 - component identification; 4 - distance of component from core 
component, in mas; 5 - position angle of component with respect to core component, in degrees; 6 - polarized flux of component, in mjy; 7 - EVPA of 
component, in degrees (nominal error of i3° from calibration, which is much greater than any model-fit errors); 8 - Degree of polarization, in per cent, 
taken from a 3x3 pixel area in degree of polarization image centred on the model-fitted position; 9 - minimum rotation measure obtained from minimum 
separation between X5 GHz and X8.4 GHzi T3.d. m"" (with nominal error of it28 rad m~ ); 10 - intrinsic EVPA, in degrees, as corrected by RMj^in. 
11 - integrated (Galactic) RM, in rad m""'^, from literature " iTaylor et al ] i2009l ). ' lOren fe Wolf j jl995l l. 12 - intrinsic minimum RM corrected for Galactic 

contribution whore RMJ^?^ = (RM„,i„ - RMGal)(l + z)^. 



thick geometrically-thin accretion disk they estimate a black 
hole mass of ~ 4 x 10^° Mq. 

Our 5 GHz image shows a core-jet structure extend- 
ing Southwards to ^ 11 mas, c onsistent with the VSOP 
Space- VLBI image at 1.6 GHz of iHirabavashi et all (|l998l ). 
which has about a factor of two worse resolution than our 
own image in the North-South direction. At 8.4 GHz the 
extended jet emission is very faint. If we cross-identify an 
intensity component in the inner jet with the innermost jet 
component at 5 GHz in order to align our images, we ob- 
tain a physically reasonable spectral index map. After ap- 
plying this same relative shift to the polarization model fits, 
both fits indicate a polarized component close to the ori- 
gin, with a residual offset between frequencies of less than 



0.20 mas; we accordingly have identified both components 
with the core polarization. A comparison of the correspond- 
ing model-fitted core EVPAs indicates a minimum RM of 
384 rad m~^ 



3.2 0636+680 

This source has been o bserved previo u sly on mas scales in 
the radio at 5 GHz bv iGurvits et aP ()l994h and more re- 
cently bv lFev fc CharlotI l|l997h with the VLBA at 2.3 and 
8.6 GHz who found it to be unresolved. Their resolution at 
8.6 GHz (1.3 X 1.2 mas) is similar to our resolution at 5 GHz 
(1.25 X 0.91 mas). It was first reported as a GPS source by 
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lO'Deal (I1990I'). The redshift of 3.17 is quoted from the cata- 
logue of IVeron-Cettv fc VerorJ (|l989l ). 

Due to the lack of extended jet emission at 8.4 GHz the 
images could not be aligned based on their optically thin jet 
emission. Polarized emission was detected in the core in our 
5 GHz map at the level of m ~l-2%; no polarized emission 
was detected at 8.4 GHz above 0.2 mjy. 



3.3 0642+449 

This extremely luminous quasar (Lg.4GHz = 6.8 x 10^* W 
Hz~^; the highest in our sample) is regularly observed by 
the MOJAV^B team with the VLBA at 15 GHz, who have 
reported su bluminal speeds o f P ~ 0.76 for an inner jet 
co mponent (iLister et al.ll2009l ). A 5-GHz global VLBI image 
bv lGurvits et al.l l|l992! T shows the jet extending almost 10 
mas to the East while |Volvach (2003) find the jet extending 
~ 5 mas to the East from EVN observations at 1.6 GHz. 

Our 5-GHz im age i s consistent with t he images of 
iGurvits et all (|l992l ) and lVolvach et al.l (|2000l ). with no ex- 
tended jet emission detected at 8.4 GHz. Due to the lack 
of extended jet emission at 8.4 GHz the images could not 
be aligned based on their optically thin jet emission. Polar- 
ized emission is detected in the core at both 5 and 8.4 GHz; 
the minimum difference between the model-fitted polariza- 
tion angles is less than 1°, nominally indicating negligible 
Faraday rotation. 



3.4 1351-018 

This radio loud quasar is the third most distant sourc e in 
our sample with a redshift of 3.707 (|Osmer et al.l 119941 1. It 
was o bserved by the V SOP Space- VLBI project at 1.6 and 5 
GHz (|Frev et al.ll2002h . It has a complex pc-scale structure 
with the jet appearing to bend through ~ 130°. The high 
resolution space- VLBI images (1.56 x 0.59 mas at 5 GHz) 
clearly resolve an inner jet component within 1 mas of the 
core with a po sition angle of ~ 120° , which was also detected 
at 8.6 GHz bv lFev fc CharlotI (|200G| ). 

Our 5-GHz image is dominated by the core emission, 
but the polarization gives a clear indication of the presence 
of an inner jet component, which is also indicated by the 
intensity model fitting. Although two polarized components 
are visible in the 8.4-GHz map, the relationship between 
these and the total-intensity structure is not entirely clear. 
The polarization angle of the polarized feature slightly north 
of the core is similar to the EVPA of the 5-GHz core, and we 
have on this basis identified this polarization with the 8.4- 
GHz core. The beam is relatively large in the North-South 
direction, and it is possible that this has contributed to the 
shift of this polarization from its true position relative to 
the intensity structure. If this identification is correct, the 
minimum RM obtained for the core is —223 rad m~^. 



complex, the intensity structures match well at the two fre- 
quencies. The components in the innermost jet lie along a 
position angle of —24°, con sistent with the higher resolution 
image of I Yang et all (|2008l ). who find an inner jet direction 
of -26° from 5-GHz VSOP data. 

Polarized components A, B, C, and D agree well be- 
tween the two frequencies, although it is not clear how these 
correspond in detail to the intensity components in the same 
regions. There is an additional polarized component Al in 
the 8.4-GHz core, which does not have a counterpart in in- 
tensity or at 5 GHz; it is difficult to be sure whether this is 
a real feature or an artefact. 

The core polarization structure at both 5 and 8.4 
GHz shows three distinct features within the core region. 
If we compare th is with the VSOP space- VLBI image of 
lYang et all (|2008l ) at 5 GHz we see that a similar type of 
structure is seen in total intensity. This indicates how the 
polarized emission can give information about the jet struc- 
ture on scales smaller than is seen in the / image alone. 



3.6 1508+572 

This is the most distant obje ct in our sample at a red- 
shift of 4.30 |Hooket^[l99i). Hence, the frequencies of 
the emitted synchrotron radiation are the highest in our 
sample at 24.6 and 44.6 GHz. The images presented here 
are the first observations to show the direction of the inner 
pc-scale jet. This quasar also has an X-ray jet extending in a 
south-westerle y direction on kpc (arcsecond) scales detecte d 
with Chandra jSiemiginowska et al.ll2003 : lYuan et al.ll200j ) . 
iGheund (|2004l ') detected a radio jet in the same region using 
the VLA at 1.4 GHz. 

Although an optically thin jet component is detected 
at both frequencies roughly 2 mas to the south of the core, 
this component proved too weak to be used to align the two 
images. Both the core and this jet component were detected 
in polarization at 5 GHz, while only the core was detected 
at 8.4 GHz. The smallest angle between the 5 and 8.4 GHz 
core EVPAs is only 1°, which we take to indicate an absence 
of appreciable Faraday rotation. 



3.7 1557+032 

This quasar is the seco nd most distant object in our sample 
with a redshift of 3.90 (|McMahon et al.lll994l '). 

There are two distinct polarized components in the core 
region where only one total intensity component is distin- 
guished; the fact that these two components are visible and 
model fit at both frequencies with similar positions and po- 
larization angles suggests that they are real. This appears to 
be a case when the polarized emission provides information 
on scales smaller than those evident in the total intensity 
image. 



3.5 1402+044 

This flat-spectrum radio quasar has a redshift of 3.208 
l|Barthel et al.lll990l ). Although the source structure is fairly 



http: / / www.physics.purdue.edu / astro/MO JAVE /index. html 



3.8 1614+051 

This quasar, at a redshift of 3.212 jBarthel et al.lll990l ') , has 
been observed bv Jev fc CharlotI (|2000l ') with the VLBA at 
2.3 and 8.6 GHz, who observed the jet to lie in position 
angle —158°. Our 8.4 GHz image shows the source clearly 
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resolved into a core and inner jet, extending in roughly this 
same position angle. 

Both our 5 and 8.4-GHz data are best fit with two com- 
ponents, whose positions agree well at the two frequencies. 
One polarized component was model-fit at 5 GHz (Table 
4); its position does not completely agree with the positions 
of either of the two / components, although it clearly cor- 
responds to jet emission. Polarization from both of the / 
components is detected at 8.4 GHz (Table 4), as well as 
another region of polarized emission between them. Com- 
ponent C in the 8.4-GHz polarization fit is weak, but its 
position agrees with that of the jet I component fit at this 
frequency, suggesting it may be real. 



3.9 2048+312 

IVeron-Cettv fc VeronI (jl993l ) found a redshift of 3.185 for 
this quasar. The apparent shift in the position of the peak 
between 5 and 8.4 GHz is an artefact of the mapping process; 
the model fits indicate a core and inner jet whose positions 
agree well at the two frequencies, as well as another weaker 
jet component further from the core detected at 5 GHz. 

This is a promising candidate source for follow-up multi- 
frequency polarization observations because the jet is well 
resolved with VLBI and has a strongly polarized core and 
jet. 



3.10 2215+020 

iDrinkwater et al] lll997) found a n emission redshift of 3.572 
for this quasar. iLobanov et al.l ((2001,) present a 1.6 GHz 
VSOP space- VLBI image of this source showing the jet ex- 
tending to almost 80 mas (~ 600 pc) with a particularly 
bright section between 45 and 60 mas. The extent of this jet 
is ~ 4 times greater than in any other pc-scale jet observed 
for quasars with z > 3. 

Our 5 GHz im age has similar resolution to the 
iLobanov etaP (|200ll ) image but is less sensitive to the ex- 
tended jet emission hence, our image has a similar but 
sparser intensity distribution. We also detect the particu- 
larly bright region, where the jet changes direction from East 
to North-East on the plane of the sky. 



4 DISCUSSION 

4.1 Brightness Temperature 

The median core brightness temperature at 5 GHz is 5.6 x 
10^^ K, while at 8.4 GHz, the median value is slightly smaller 
at 3.2 X 10^^ K. This can also be seen from a histogram of the 
core brightness temperature (Figure |3} where the 8.4 GHz 
values populate a larger majority of the lower bins. This may 
be as a result of the 8.4 GHz data probing regions of the 
jet where the physical conditions are intrinsically different, 
leading to lower observed brightn ess temperatures, similar 
to the results of iLee et all l|2008l ') at 86 GHz. However, it 
is difficult to separate this effect from any bias due to the 
resolution difference between 5 and 8.4 GHz. Furthermore, 
due to the relatively small difference between the median 
values, a larger sample size would be required to determine 



whether this is a real effect or just scatter in the data for 
the small number of sources observed here. 

The maximum core brightness temperature at 5 GHz 
is 4.2 X 10^^ K found in 0642+449 and the minimum value 
is 3.7 X 10^" K found in 2048+312. At 8.4 GHz, the max- 
imum and minimum values are 7.3 x 10^^ K (0642+449) 
and 7.4 x 10^" K (0636+680). Assuming that the intrinsic 
brightness temperature does not exceed the equipartition 
upper limit of Teq ~ 10" K (|Readheadlll"993 ). we can con- 
sider the observed core brightness temperatures in excess 
of this value to be the result of Doppler boosting of the 
a pproaching jet emission. Us ing the equipartition jet model 
of lBlandford fc Konigll (|l979D for the unresolved core region, 
we can estimate the Doppler factor (5) required to match the 
observed brightness temperatures (Tb,obs ~ 3 x IG^^S^^^). 
In most cases the Doppler factors required are modest, with 
values ranging from 1 to 5, except for 0642+449 which re- 
quires a Doppler factor of 23 at 5 GHz and 46 at 8.4 GHz. 

From the observed brightness temperatures of the in- 
dividual jets components, we can invest igate the a ssump - 
tion of adiaba t ic ex pansio n following Marscher ([l99y) , 
ILobanov et al.l |2003) and iPushkarev et al. I l|2008l ). who 
model the individual jet components as independent rela- 
tivistic shocks with adiabatic energy losses dominating the 
radio emission. Note that this description of the outer jet dif- 
fers from the N oc , B oc case describing the compact 
inner jet region, which is not adiabatic. With a power-law 
energy distribution of N{E) cx E^"'^ and a magnetic field 
described by B oc d^", where d is the transverse size of the 
jet and a = 1 or 2 corresponds to a transverse or longitudinal 
magnetic field orientation, we obtain 



7b, 



2~b core 



encore 
d\ct 



a + l-ci(a + 4/3) 



(3) 



which holds for a constant or weakly varying Doppler fac- 
tor along the jet. Hence, we can compare the expected jet 
brightness temperature (Tb.jet) with our observed values 
by using the observed core brightness temperature (Tb^coro) 
along with the measured size of the core and jet components 
(dcorc and djct) for a particular jet spectral index. We apply 
this model to the sources with extended jet structures to see 
if this model is an accurate approximation of the jet emis- 
sion and also as a diagnostic of regions where the physical 
properties along the jet may change. 

Figure 13] shows the comparison of this model (solid line) 
with the observed brightness temperatures (dashed line) 
along the jet of 0014+813 at 5 GHz. From the spectral in- 
dex distribution obtained between 5 and 8.4 GHz, we adopt 
a jet spectral index of a = —1.0 and assume a transverse 
magnetic field orientation (a = 1) from inspection of the jet 
EVPAs (Fig. 2). The two jet components furthest from the 
core agree relatively well with the model, however, the two 
jet components closest to the core are substantially weaker 
than expected. Given the high RM for this source, it is also 
possible that the intrinsic magnetic field orientation is lon- 
gitudinal; using a value of a = 2 yields model brightness 
temperatures for these components that are within our mea- 
surement errors. Unfortunately, the extended jet emission at 
8.4 GHz is not detected so we cannot constrain the jet RM 
and we also cannot rule out a change in the Doppler fac- 
tor along the jet due to either an intrinsic change in the 
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jet dynamics or a change in the jet direction. Clearly there 
are not enough observational constraints for this source to 
convincingly test the applicability of this model. 

The observed brightness temperatures in the jet of 
1402+044 at both 5 and 8.4 GHz (Figure[5| are in reasonable 
agreement with the model predictions for a = —0.5, which 
is consistent with the overall jet spectral index distribution. 
Since the EVPA orientation implies a transverse magnetic 
field (Figure [2]), we use B cx . The jet components at 
1 mas and 9 mas from the core both have observed bright- 
ness temperatures slightly higher than would be expected 
for radio emission dominated by adiabatic losses. Using flat- 
ter spectral index values of a = 0.0 and a = —0.3, respec- 
tively, for these two components brings the model brightness 
temperatures within the measurement errors. Hence, it is 
likely that these components are subject to strong reaccel- 
eration mechanism that temporarily overcomes the energy 
losses due to the expansion of the jet. 

Applying the analysis to the jet of 2215-)-020 at 5 GHz 
(Figure |6} , we see that the model matches the observed val- 
ues very well for the two inner jet components using our mea- 
sured spectral index of a = —0.7 and a transverse magnetic 
field (a = 1). However, the extended components brightness 
temperatures at 50-60 mas from the core are higher than 
expected by the adiabatic expansion model. This suggests 
different physical conditions in this region of the jet creating 
a less steep spectral index (using a = —0.5 brings the model 
in line with the observed values). This far from the core it 
is likely that the Doppler factor has changed with either a 
change in the jet viewing angle or speed, possibly from inter- 
action with the external medium. This is consistent with the 
observed strong brightening of the jet along with its change 
in direction from East to North-E ast, shown in Figured Our 
results are very similar to those of lLobanov et al] (|200ll ) who 
employed the same model for this source at 1.6 GHz. 

4.2 Jet Structure and Environment 

The added value of polarization observations in providing 
information on the compact inner jet structure is clear from 
the images of 1402-1-044 (Figure [2]), 1557-1-032 (Figure [g]) 
and 1614-1-051 (Figure [2|, where polarized components are 
resolved on scales smaller than obtained from the total inten- 
sity image alone. For 1402-1-044, a higher resolution VSOP 
image at 5 GHz (jYang et al.l l2008l ) resolved three compo- 
nents in the compact inner jet region consistent with what 
we find from the polarization structure of our lower reso- 
lution images where only one total intensity component is 
visible. In the case of 15574-032 and 1614-1-051, the iden- 
tification of multiple polarized components in the core re- 
gion, where again only one total intensity component is 
directly observed, suggests that observations at higher fre- 
quencies and/or long er baselines p rovided by, for example, 
the VSOP-2 mission (lTsuboi|[2009l ) are likely to resolve the 
total intensity structure of the core region. Hence, the po- 
larization information is crucial in identifying the true ra- 
dio core and helping to determine whether it corresponds 
to the T = 1 frequency-dep endent surface first proposed by 
iBlandford &: K onigl (1979) or a stationary feature such as 
a conical shock (Cawthornc 2006; Marscher 2009 ) for a par- 
ticular source. 

For the sources with extended jets, substantial jet bends 



are observed. This is not a surprise since as long as the ob- 
servations imply relativistic jet motion close to the line of 
sight, small changes in the intrinsic jet direction will be 
strongly amplified by projection effects. For example, an 
observed right-angle bend could correspond to an intrinsic 
bend of only a few degrees (e.g., ICohen et al.ll2007h . In the 
case of 1351-018, the jet appears to bend through ~ 130° 
(Figure [2]) from a South-Easterly inner jet dir ection to a 
North- Westerly extended jet direction (see also iFrev et al] 
12002"). The polarization distribution in the core region also 
supports an South-Easterly inner jet direction. This may be 
due to a helical jet motion along the line-of-si ght as pro- 
poscd, for example, for the jet of 1156-1-295 by iHong et al.l 
(2004) or small intrinsic bends due to shocks or interactions 
with the external medium amplified b y projection effects i n 
a highly beamed relativistic jet (e.g., iHoman et al.ll2002l ). 
The jet bend at ~ 400 pc from the core of 2215-1-020 at 
5 GHz may be due to a shock that also causes particle reac- 
celeration, increasing the observed jet emission. 

The core degree of polarization found in these sources 
is generally in the range 1-3%. Given that the emitted fre- 
quencies are in the range 20-45 GHz, we can compare these 
values to the core degrees of polari zation observed in low- 
redsh ift AGN jets at 22 and 43 GHz (|0'Sullivan fc Gabuzdal 
120091 ). In general, the core degree of polarization is higher 
in the low-redshift objects with values typically in the range 
3-7%. This could be due to intrinsically more ordered jet 
magnetic field structures at lower r edshifts or more likely, 
selection effects, since the AGN in lO'Sullivan fc Gabuzdal 
(|2009l ) were selected on the basis of known rich polariza- 
tion structure while the high-redshift quasars had unknown 
polarization properties on VLBI scales. Another likely pos- 
sibility is that the effect of beam depolarization (different 
polarized regions adding incoherently within the observing 
beam) is significantly reducing the observed polarization at 
5 and 8.4 GHz due to our much larger observing beam com- 
pared with the VLBA at 22 and 43 GHz (factor of -^3 better 
resolution even considering the longer baselines in the global 
VLBI observations). 

Due to the extremely high redshifts, an observed RM 
of 50 rad m~^ at z ~ 3.5 (average redshift for the sources 
in our sample) is equivalent to a nearby source with an RM 
of 1000 rad m~^. The mean core |RMIJ5[n| for our sample is 
5580±3390 rad m~^. Naively comparing ou r results with the 
8 -15 GHz s ample of 40 AGN reported in IZavala fc Tavlod 
(120031 . |2004| ). which has a median intrinsic core RM of ~ 400 
rad m~^ (and a median redshift of 0.7), suggests that the 
VLBI core RMs are higher in the early universe. However, 
RMs of the order of 10* rad m~'^ have been measured in 
low to me dium redshift sources at similar em i tted f requen- 
cies (e.g., lAttridge et al. I l2005l : lAlgaba et af] |2010D : while 
some of the largest jet RM estimates, of order 10^ rad m~'^, 
have come through observation from 4 3-300 GHz of source s 
with redshifts ranging from - 0.1-2 (jjorstad et al.ll2007l ). 
Furtherm ore, 15-43 GHz RM m eas urements of BL Lac ob- 
j ects f rom lGabuzda et all (|2006l ) and lO'Sullivan fc Gabuzdal 
(|2009D have a median intrinsic RM of ~ 3000 rad m ^ and 
a median redshift of 0.34. Hence, from the small sample of 
minimum RM values presented in this paper it is not clear 
whether these sources are located in intrinsically denser en- 
vironments or if the large RMs are simply due to the fact 
that the higher emitted frequencies are probing further up- 
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stream in the jet where it is expected that the electron den- 
sity is higher and/or the magnetic field is stronger. 

Our results from 5 to 8.4 GHz generally show flat to 
slightly steep core spectral indices with a median value of 
a = —0.3. However, there is a large range of values going 
from 0.7 for 0642+449 to -1.0 in the case of 1614+051. In 
the majority of cases, we find steep jet spectral indices rang- 
ing from a ~ —0.5 to —1.5 with a median value of a = —1.0. 
High-redshift objects are often searched for through their 
steep spectral index due to the relatively higher emitted 
frequencies for a particular observing frequency. The z~a 
correlation, as it is known , (i.e., steeper a at hig her z) has 
been very successful (e.g.. [Broderick et al. I l2007l. and refer- 
ences therein) in finding high-redshift radio galaxies. While 
steeper spectral index values are expe cted simply fr o m the 
higher rest frame emitted frequencies, iKlamer et al.l (l2006l ) 
show that this cannot completely explain the correlation. 
The possible physical explanation they present is that the 
sources with the steepest spectral index values are located 
in dense environments where the radio source is pressure 
confined and hence, loses its energy more slowly. This effect 
might also apply to the dense nuclear environments of our 
sample of quasar jets. To test this we have plotted the core 
spectral index versus |RMJJJf^| in Figure [7] The data suggest 
that this is a useful avenue of investigation with a Spear- 
man Rank test giving a correlation coefficient of -0.95, but 
clearly more data points are needed to determine whether 
there truly is a correlation or not. 

Further observations of these sources with much bet- 
ter frequency coverage is required to properly constrain the 
spectral and KM distributions and to make detailed com- 
parisons with low redshift sources to further investigate the 
quasar environment through cosmic time. 



5 CONCLUSION 

We have successfully observed and imaged 10 high-redshift 
quasars in full polarization at both 5 and 8.4 GHz using 
global VLBI. Model fitting two-dimensional Gaussian com- 
ponents to the total intensity data enabled estimation of the 
brightness temperature in the core and out along the jet. The 
observed high core brightness temperatures are consistent 
with modestly Doppler-boosted emission from a relativistic 
jet orientated relatively close to our line-of-sight. This can 
also explain the dramatic jet bends observed for some of our 
sources. 

The core degrees of polarization are somewhat lower 
than observed from nearby AGN jets at similar emitted fre- 
quencies. However, beam depolarization is likely to have 
a stronger effect on these sources compared to the higher 
resolution observations of nearby sources. Model- fitting the 
polarization data enabled estimation of the minimum RM 
for each component since obtaining the exact value of the 
RM requires observations at three or more frequencies. For 
sources in which we were able to remove the integrated (fore- 
ground) RM, we calculated the minimum intrinsic RM and 
found a mean core |RMUJf^| of 5580 rad m~^ for four quasars 
with values ranging from —1249 rad m~^ to 9110 rad m~^. 
We found relatively steep core and jet spectral index val- 
ues with a median core spectral index of —0.3 and a me- 
dian jet spectral index of —1.0. We note that the expecta- 



tion of denser environments at higher redshifts leading to 
larger RMs can also lead t o steeper spectral in dices through 
strong pressure gradients (|Klamer et al.|[2006l ) and that this 
hypothesis is not inconsistent with our results. 

Several of the sources presented in this paper are inter- 
esting candidates for follow-up multi-frequency observations 
to obtain more accurate spectral and RM information in or- 
der to make more conclusive statements about the environ- 
ment of quasar jets in the early universe and to determine 
whether or not it is significantly different to similar lower 
redshift objects. 
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V vs U for 2048+312 at 5 GHz 
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Figure 1. Example plot of the uv coverage obtained for a representative source, 2048+312 at 5 GHz. 
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Figure 2. Total intensity (7) and polarization (p) maps for all sources at 5 GHz (left) and 8.4 GHz (right). The image axes show the 
relative offset from the image centre in milliarcseconds. The beamsize is indicated by an inset in each image (see Table 1 for details). 
The model- fitted total intensity components are indicated by the "+" symbols while the position and angle of the polarized component 
model-fits are indicated by the boldface line. In all cases the contour levels increase in factors of 2 from the indicated total intensity 
cutoff (see Table 5). The length of the p sticks represent the relative polarized intensity. 
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Figure 2. continued. 
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Figure 2. continued. 
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Figure 3. Histogram of core brightness temperature, in units 
of at 5 GHz (grey) and 8.4 GHz (white) for 9 sources 

(0642+449 excluded). 
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Figure 4. Brightness temperature (Tf,) for components along the 
jet of 0014+813 at 5 GHz. The points connected by the dashed 
line represent the observed brightness temperatures while the 
points connected by the solid line represent the model values. 



1402+044: 5 GHz 
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(a) 5 GHz data. Dashed line: observed; Solid line: model predietion 
with a — —0.5. 



1402+044: 8 GHz 
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(b) 8.4 GHz data. Dashed line: observed: Solid line: model prediction 
with a — —0.5. 

Figure 5. Brightness temperature (Tj) for components along the 
jet of 1402+044 for both 5 and 8.4 GHz. The points connected by 
the dashed line represent the observed brightness temperatures 
while the points connected by the solid line represent the model 
values. 



2215+312: 5 GHz 
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Figure 6. Brightness temperature (Tj,) for components along the 
jet of 2215+312 at 5 GHz. The points connected by the dashed 
line represent the observed brightness temperatures while the 
points connected by the solid line represent the model values. 
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Figure 7. Plot of the intrinsic minimum RM (RMJ^[^) versus 
the core spectral index between 5 and 8.4 GHz for each source. 
A Spearman Rank test gives a correlation coefficient of -0.95. 
However, more data points are required to determine whether 
there is a correlation or not. 
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